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. . INTRODUCTION 


The operation of a rocket engine at various pressure regions in 
the combustion chamber results in an uneven side-load in the rocket 
nozzle before the engine reaches the full operation condition. One 
solution proposed for preventing the uneven side-load due to the separ- 
ation of flow, is to use the "trip wires" in the nozzle to create an even 
forced flow separation and to reduce side-load during the start transient. 
The purpose of the present study is to investigate the characteristics of 
flow profile in the rocket nozzle during the start transient, and to 
examine the possibilities of reducing the side-load thrust by sticking 
"trip wires" in the nozzle. 

To simplify the geometry of flow configuration around the trip 
wires, we have assumed that the flow is passing through square steps 
instead of round wires. Since a purely analytic solution is not availa- 
ble at this moment, we have proposed a series of semi-empirical solutions 
based on the existing literature. In Chapter II, the possible mathe- 
matical formulations of the flow parameters, which describe the flow 
patterns in front of and after the square step, are presented. In 
Chapter III, the boundary layer thickness and the values of accompanied 
flow parameters through whole regions of specified nozzle are calculated 
numerically. In this calculation, tb interaction between boundary layer 
and shock wave has been ignored. The value of boundary layer thickness 
at the appropriate location of the nozzle is a significant factor for the 
delsgn of step height in the experimental investigations. 

Experimencal setup and the model of nozzle used in our experiment, 
with appropriate steps mounted on the wall of the nozzle, are described 
in Chapter IV. The experimental results and the discussions are presented 
in Chapter V. 


II. POSSIBLE MATHEMATICAL MODELS 


As stated earlier, the nozzle trip wire configuration was proposed 
to be a potential method for reducing start- transient side loads. Theo- 
retical analysis shows that the side loads are apparently a result of the 
unsteady, non-uniform flow separation from the wall of the nozzle. Thus, 
an installation of a device to force a uniform separation should result 
in a reduction of the magnitude of the side loads. However, a single 
trip wire would provide the mechanism to control separation at one specific 
point in the nozzle. As the amount of separation realized depends on the 
chamber pressure level, it is apparent that more than one trip wire would 
be required. Therefore, a system of multi-trip wires was chosen. 

The basic concept of multi-trip wires is the ability of the flow 
to snap from one wire to the downstream wire and yet provide a shock of 
sufficient strength to maintain uniform separation providing the mechanisms 
for reducing the transient side loads. The problem is how to choose trip 
wire size and spacing required for each trip wire. Theoretically, these 
problems are concerned with the study of boundary layer-shock wave inter- 
action in the flow separation associated with recompression in a super- 
sonic nozzle. 

In general, boundary layer approach, using continuity, momentum, 
and energy equations, is used extensively by investigators to solve 
classical flow separation problems. 

5c + 57 + £ (t V • 0 (2 ' 1) 



for compressible two dimensional boundary layer. Here p, u, P, T, and C p 
denote fluid density, velocity, pressure, temperature and constant 
pressure specific heat, respectively. Subscripts x and y Imply component 
along x and y directions. Coefficients K and y ere thermal conductivity 
and viscosity. A complete culation of the boundary layer for a given 
body with the aid of Equations (2-1) end (2-3) is, in many cases, so 
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cumbersome and time consuming that it cannot be carried out in practice. 
It is, therefore, desirable to possess at least approximate methods of 
solution, to be applied in cases when an exact solution of the boundary 
layer equations cannot be obtained with a resonable amount of work, 
even if their accuracy is only limited. Thus, a mean value is obtained 
from the momentum equation which is derived from Equation (2-2) by inte^ 
gration over the boundary layer thickness. The equation is known as the 
momentum integral equation of boundary layer theory, or as Von barman’s 
momentum integral which is 


Cb 


— * -J- (U 2 6**) + 5* U 
p dx » » dx 

for two dimensional Incompressible boundary layers. Here 


(2-4) 


/ 3u\ 

t - p l “j” ) * shearing stress at the wall 

' 'y-0 


! 0 for y ■ 0 (at wall) 

U y * 5 (outside the boundary layer) 


5 * 



(displacement thickness) 


6 ** « 



-■7- J dy (momentum thickness) 


To deal with the interaction between an internal dissipative flow and an 
external nearly isentropic stream such as shock wave interaction with 
boundary layer, Crocco and Lees 1 generalised the concept of the Von 
Kerman momentum Integral for the dissipative flow region. Mixing between 
the external and Internal streams Is admitted as a fundamental physical 
process, which will furnish the mechanism for pressure rise. This mixing 
theory is applicable not only for separated and reattached flows, but 
also for wake flows. One of the most significant discoveries by Crosco 
and LeeS 1 Is that In the separation flow the pressure gradient along the 
surface may reach a maximum at separation and may drop off steeply down- 
stream, and that in reattaching of wake flows the pressure gradient is 
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negligible some distance upstream of the reattachment point and in- 
creases rapidly as this point is approached. In flow separation two 
governing parameters k and f are defined on the basis of the boundary 
layer as (Crocco and Lees 1 ) 

k * 

and 

f * 

where the subscript 1 denotes the incompressible flow in the Stewartson 
transformation (Chang 2 ). 

For laminar and turbulent flows, the separation is characterized 
by decreasing k and by increasing f. At the separation point, the slope 
df/dk is infinite. This physical picture is essential for separated flow 
solutions since it gives a clue to the quantitative prediction of what 
happens after separation. No reliable theoretical solution 
has been established in this region. 

In the practical application, it is assumed that the wire is 
attached to a flat wall and that the flow is separated on both the up- 
stream and downstream sides of the wire (see Figure 1). Since the wires 

are small relative to the nozzle diameter, the flow in the region of the 
wire can be assumed to be planar. Of course, the flow from the throat 
of the nozzle to the region of first wire is assumed to follow the theory 

of isentropic expansion. A single oblique shock is taken to occur at 

the beginning of the front separated region and a second at the downstream 
edge of the separated base region. Prandtl-Meyer expansion fans were 
considered to occur at the upper two corners of the wire. The height of 
the wire is chosen to be comparable to the thickness of the boundary 
layer. Calculation of the boundary layer thickness In the present study 
will be shown In the next chapter. 

To make the mathematical model tractable, the simplified formulations 
shown by Chang 2 have been used. 

A. Front Step 

The pressure coefficient on the front side of the wire is 



6 - 6 * - 6 
6 - 6 * 
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(2-7) 





Flow Conf igur otion Around The Step 









where the subscript denotes regions of flow under considerati 
figure 2 ). 

The relationship among C^, shock angle resulting fro. 
region 8 * and angle of boundary of separated region 0 ^, is 


( ^ (see 
.separate 


1 c cot 6 

tan 0 , * • ( 2 - 8 ) 

2 "p 

Thus, the value of the angle 0^ can be obtained if the shock, angle 8 can 
be measured from the experiment. The position of the boundary layer 
separation ahead of the step can be obtained from the following relation 


= h cot 0 £ 

where h is the height of step. 


(2-9) 


B. Back Step 

The flow across the regions and can be calculated based on 
Prandtl-Meyer expansion 


V 2 ( V + e 2 = V 3 (M 3 } + 9 3 


(2-10) 


where 6 ^ * 0^ and 0^ * 0. Thus, Equation (2-10) becomes 


V 3 (M 3 } “ V 2 ( V + 6 f * 


( 2 - 11 ) 


This provides the Mach number M^. Furthermore, from momentum equation we 


can calculate the pressure as follows: 


P 3- P 2 


1 + YM‘ 

1 + y m* 


( 2 - 12 ) 


Again flow across the regions and shall be governed by Prandtl- 
Meyer expansion. If we use the numerical results computed by Korst 3 in 
which the results show a good agreement among the experimental data made 
by Chapman, et .>1.* for the pressure ratio P^/P^ against the 
approaching Mac', number (M^) (See Figure 3), we can calculate the Mach 
number based on the following relation 
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p* O+YMO - 1 
4 


(2-13) 


Then, Che angle of boundary of separated region 0^ can be obtained from 
Prandtl-Meyer expansion 


\ = V 4 * M 4^ ” V 3 (M 3 ) 


(2-14) 


Therefore, the position of boundary layer separation behind the step 
can be calculated as follows: 


•= h cot 6^. 


(2-15) 


Through experiment curve fitting, Herbert and Herd* suggested 


4 2 

— - 1.15 + 0.05 M, . 

p. 3 


Combining with P,/P,» as in Figure 3, a plot of S./h against M, is 

J *4 D J 

shown in Figure 4. All notations are defined as in Figure 2. Figure 4 
shows that the length of the reattachment decreases with increasing 
Mach number. 
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III. COMPUTATION OF BOUNDARY LAYER THICKNESS 


The computation of boundary layer, in general, shall follow the 
theory of boundary layer developed by Prandtl. However, Von Karman 
integral momentum equation is more suitable for practical application 
to engineering problems. L. this section, we use the integral momentum 
and energy equations in axisymmetric form for compressible boundary layer 
flow. The governing equations are as follows: 


d5** _ f , , /dr \ 2 

dx 2 \dx/ 


+ i r- 

r dx 


for integral momentum equation, and 


1+ dU . d(p U ) 

A ** 6** ® , 1 00 00 

U dx p U dx 

oo oo 00 


(3-1) 


JA H - H 

d4> _ r aw w 

dx H H - H 
o 


J_ d(p - U . ) j. dr 1 d( VV 

U dx r dx H - H dx 

to oo O W 


for integral energy equation. Here, the displacement thickness 5*, 
momentum thickness 6** and energy thickness $ for the compressible flow 
are identified as follows: 


■ \ ( i "&) 
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The skin friction coefficient is defined as 

2t 


C f = 


P u 

00 oo 


which has a form of Blasius relation 


c f = 


0.025 


( Re 6 „) 0 * 25 


where the Reynolds number in terms of momentum thickness is 


Re x** 


p U (5 
00 00 


** 


The Stanton number is defined as 


C H " p 0 (H -H ) 
oo » aw w 

which can be shown in the following form 


c f (I V _ (a 


C H = 


r ** 


C f (Re 4 ) 1/2 6 

1-5 - ..-I- i - P + Jin — 

2 r 5P +1 


Velocity and enthalpy profiles across the boundary layer are assumed to 
follow the relationships: 


for 


h - H 
o w 

H - H 
o v 


■(f) 


for y < 6 


• 1 for y > 6 
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h -H 
o w 

H -H 
o w 


= 1 for y > 6 


where o and are momentum and thermal boundary layer thicknesses, 
respectively. The definition of enthalpy is 


H = I C d T 

I P 


h = H+ -V 
o 2 


w 

W 


C dT 
P 


The adiabatic wall enthalpy H is defined as 

3W 


H „ + ( P r > l/3 T 


H » + T 


The density p within the boundary layer is obtained from the perfect 
gas equation, assuming the pressure is constant across the boundary 
layer 


P m T ' 

The computer calculation is carried out through the programming made 
by Omori et al. *. The results are shown in Figure 8. 


■ Ml 




IV. APPARATUS AND MODELS 


A. Apparatus 

The set up of the experiment is shown in Fig. 5 and Fig. 6. A 
schematic diagram of the flow system is shown in Fig. 7. Flow pictures 
were taken by using a differential interferometer which is also shown 
in Fig. 7. The test section was 1.8" x 2.1". The compressed air was 
stored in a 300 cubic feet tank which allowed the reservoir pressure 
to be maintained at 600 psig during the entire test. 

B. Experimental Model 

The shape of the model is shown in Fig. 8. It is 2.86" long and 
0.8" wide. The upper part of the nozzle was fixed. In order to change 
step heights easily for different tests, the lower part of the model was 
made removable. Slots of depth of 0.025" were cut in the model for the 
step to be mounted. The location of the slots are 0.75", 1.055", and 
1.358" from the end of the nozzle exit, respectively. 

The heights of the step were made 0.027, 0.033, 0.035, 0.039, 
0.045, 0.051, and 0.057 inches which gave the step heights 0.002, 0.006, 
0.008, 0.014, 0.020, 0.026 and 0.032 inches when plugged into the slot 
of the model. Thus, we were able to use different step heights between 
0.002 and 0.032 inches, which covered the range of both lower and higher 
than the boundary layer thickness In the test sections. 

In the later test of two or three step experiments, models were 
cut two or three slots at proper position chosen. All possible bleed 
passages were sealed with sealant. A small air bleeding was used to 
prevent formation of fog on the window of the test section. 
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V. TEST PROGRAM 


In order to simulate the start transient condition of rocket 
engine in the present test program, the chamber pressure was increased 
from zero psig to normal operation pressure up to about 300 psig which 
allowed expansion nozzle exit Mach numbers up to ~ 3.5. The test pro- 
cedure is to introduce the high pressure gas into chamber gradually 
and increase the gas flow rate paslng through the expansion nozzle. 
During the test, we have observed how the shock was formed and passed 
over the steps, and how it propagated out of the nozzle when both gas 
pressure and flow rate were increased. The flow parameters along the 
nozzle, such as temperature, pressure, density, velocity, Mach number 
and boundary layer thickness, are shown in Fig. 9. Photographs were 
taken in order to determine the natural separation points at various 
pressures. These are needed for the calculation of the ratio of 
separation pressure to ambient pressure. 

Single step experiments were carried out by inserting various 
step heights into the slots at three different positions as previously 
described. For each test case, reservoir pressure at the moment the 
shock attached to the step and at the moment the flow reattached after 
flow separation, were recorded. The purpose of these experiments was 
to determine how much increase in reservoir pressure is needed for the 
flow to reattach after being force-separated by the existence of the 
step. Tests included both change of step heights and step positions. 
Step heights used were 0.002, 0.008, 0.014, 0.020, 0.026, and 0.032 
inches. Steps were placed at 0.75, 1.055, or 1.358 inches from the 
end of the model. 

To study the change of reservoir pressures needed between flow 
separation and reattachment for multi-step conditions, two-step and 
three-step cases were also tested. In two-step cases, steps were placed 
at the positions 0.75 and 1.358 inches from the end of the model, an.', 
step heights used were 0.002, 0.006, 0.010, 0.020, and 0.032 inches. 

In three-step cases, the positions of the steps were 0.75, 1.055, and 
1.358 Inches from the end of the model and step height chosen was 0.032 
Inches. Again, the reservoir pressures were recorded at the time the 



flow separated at the step and when the flow reattached after the shock 
moved downstream. 

The next point of natural separation after the reattachment of 
flow on each step was determined photographically. 
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VI. RESULTS AND DISCUSSIONS 


Points of natural separation versus chamber pressures are shown 
in Fig. 10. The relation appears t:o be linear which means that the in- 
crease reservoir pressure results in pushing the point of natural 
separation further out. It was also found that the ratios of pressure 
at the point of separation and the atmospheric pressure are ranging 
from 0.5 to 0.41 for the range of reservoir pressure from 70 to 220 psig. 

In Table I, reservoir pressures, at the time the shock attaches 
to the step and the flow reattaches after the flow separation due to 
the step, are listed for different step heights and three different 
positions (0.75", 1.055" and 1.35" from the exit of nozzle). Local 
static pressures before the shock and after flow reattachment corresponding 
to Pqi's and ate also listed. The static pressure difference 

needed for the shock to jump over the step was found statistically to 
be about 2.5 psl. Figure 11 shows the ratio of local pressure just 
before flow separation (P^) to that of local pressure just after flow 
reattnchment (Pj) as a function of step height. Figure 11 also shows 
that this ratio Increases drastically at h/6 • 1. Table II shows the 
ratios of local pressure tc ambient pressure at the point of separation 
when the flow reattached. It is Interesting to point out that the 
pressure ratio Is about 0.9, which is much higher than the case without 
the step. From the photographic measurement, the distance from the 
flow reattachaent to the step to the next point of natural separation 
is 10 tinea the height of the step. 









Static pressure before shock as the shock attaches to the step. 

Static pressure before shock as the flow reattaches after 
separation shock moves downstream. 

Step height 

Boundary layer thickness. 


i 


1.8 



Fig. 11. P 2^ P 1 Versus Ste P Height (h/rt). 
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Step Close Step Step Close 

to Throat at Middle to Exit 


TABLE I. PRESSURE DATA FOR SINGLE STEP CASES 



P 02 

psia 

P 1 

psia 

290 

5.434 

264 

4.862 

280 

5.005 

255 

5.148 

215 

4.974 

236 

5.343 

220 

5.527 

218 

5.896 

220 

6.080 

219 

6.411 

170 

6.348 

163 

6.295 

167 

6.295 

175 

6.295 

170 

6.506 

171 

6.452 



8.294 


7.550 


8.008 


7.293 


7.922 


8.696 


8.107 


8.033 


8.107 


8.070 


8.993 


8.622 


8.834 


9.257 


8.993 


9.045 


2.86 


2.688 


3 


2.145 


2.948 


3.353 


2.580 


2.137 


2.027 


1.659 


2.645 


2.327 


2.539 


2.962 


2.487 


2.592 


Pq^: Reservoir pressure when the shock attaches to the step. 

P- 2 : Reservoir pressure when the flow reattaches after separation by 
the step. 

P. , P,: Local static pressure corresponding to P m and P ,, 

respectively. 0A ,z 

Step close to the exit: 0.75" from the exit. 

Step at middle: 1.055" from the exit. 

Step close to throat: 1.358" from the exit. 
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TABLE II. NATURAL SEPARATION PRESSURE AFTER 
REATTACHMENT OF FLOW 




p on 


p on 

ep Height 

( inchti) 

(psia) 

P a 

P a 

0.032 

136.9 

.87 

9.31 

0.026 

131.6 

.84 

8.95 

0.020 

136.6 

.89 

9.29 

0.014 

145.0 

.908 

9.86 

0.008 

142.6 

.906 

9.7 

0.002 

143.9 

.92 

9.79 


where 

P * calculated total pressure for the flow region after the 
shock reattachment. 

P = calculated local pressure at the location where natural 
separation occurred. 

= ambient atmospheric pressure. 


Tables III and IV show the variation of chamber pressures 
during the experiment for the cases of two-steps and three-steps, 
respectively. Figures 12, 13, 14, and 15 show the typical flow 
pictures of natural separation for the cases of single-step, two-steps, 
and three-steps, respectively. The photographic pictures were taken 
by using a differential interferometer. 

On the average, it is noted that within the range of Mach numbers 
tested, the distance from the point of flow separation in front of the 
step to the step is about 3 to 5 times that of the step height; and 
the distance from the step to the point of reattachment downstream is 
about twice the step height. 



I 

TABLE III. PRESSURE DATA FOR TWO-STEP CASES 


i 

i 


% 

I 


Step 

Heights 

Step 

No. 

P 01 

psia 

P 02 

psia 

P 1 



P 02 

p oi 

32 

1 

120 

170 

6.348 

9.257 

2.909 

1.46 

1000” 

2 

230 

250 

34.13 

37.09 

2.97 

1.09 

20 

■ 

120 

165 

6.348 

8.728 

2.38 

1.38 

1000" 

2 

1 

225 

245 

33.38 

36.35 

2.96 

1.09 

10 


120 

163 

6.348 

8.622 

2.274 

1.36 

1000" 

2 

220 

240 

33. *6 

34.32 

2.86 

1.09 

6 

■I 

125 

170 

6.612 

8.943 

2.381 

1.36 

1000" 

2 

200 

220 

30.36 

33.8 

2.82 

1.10 


Steps at 0.75" and 1.358" from the Exit. 
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TABLE IV. PRESSURE DATA FOR THREE-STEP CASES 


Step 

No. 

P 01 

psia 

P 02 

psia 

P 1 

psia 

P 2 

psia 

ap-p 2 -Pi 

psi 

1 

113 

195 

5.977 

10.31 

4.333 

2 

195 

229 

14.85 

18.27 

3.42 

3 

229 

295 

16.02 

20.64 

4.62 


Three Steps: 

Steps at 0.75", 1.055", and 1.358" from the exit. 
Step height 32/1000". 
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